NATIONAL ADVISORY COMMITTEE FOB AERONAUT ICS 



RECTANGULAR SHELL PLATING UiiDER 
UNIFORMLY DISTRIBUTED HYDROSTATIC PRESSURE 
By M, Neubert and A. Sommer 

Luf t f ah r t f o r s chun g 
Vol. 1?, No. 7, July 20, 1940 
Verlag von R. Oldenbourg, Munchen und Berlin 




No. 965 



Glenn L Martin College 
^Engineering and Aero- 
nautical Sciences 
Library 



Wash ingt on 
December 1940 




Ml 



NATIONAL ADVISORY COMMITTEE FOR AERONAUT ICS 1 



TECHNICAL MEMORANDUM NO.. 965 



RECTANGULAR SHELL PLATING UNDER 
UNIFORMLY DISTRIBUTED HYDROSTATIC PRESSURE* 
By M . Neubert and A. Sommer 

SUMMARY 

A check of the calculation methods used "by Foppl and 
Hencky for investigating the reliability of shell plating 
under hydrostatic pressure has proved that the formulas 
yield practical results within the elastic range of the 
material. Foppl's approximate calculation leaves one on 
the safe side. It further was found on the basis of the 
marked ductility of the shell plating under tensile stress 
that the strength is from 50 to 100 percent higher in the 
elastic range than expected by either method, 

A, INTRODUCTION 

Knowledge of the stresses in rectangular shell plating 
under hydrostatic pressure. is important to the airplane de- 
signer since he frequently has to deal with such structural- 
parts. This applies above all to hulls and floats which at 
take-off and landing or during handling in rough sea are 
subject to enormous water pressure by the waves. But then 
it also includes gasoline tanks and gasoline chambers which 
undergo considerable interned pressure because of high ac- 
celerations. Concerning the calculation there are the re- 
ports by Foppl (reference l) and Hencky (reference 2) , the 
practicability and reliability of which are checked against 
various experimental results. 

B. CALCULATION OF RECTANGULAR SHELL PLATING UNDER 
HYDROSTATIC PRESSURE 

Both authors proceed from the assumption that the 



* "Rechteckige Blechhaut unter gleichmfissig verteiltem Fliis- 
sigkeitsdruck. » Luf t f ahr t f o r s chung , vol. 17, no. 7, July 20, 
1940, pp. 207-210. 
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skin possesses no flexural stiffness (J = 0) and is not 
strained in the relaxed state. On these premises P8ppl 
has set up differential equations for thin plates with 
great deflections. His approximate solution published 
in 1920 contained, however, only the calculation for the 
square skin. A year later Hencky gave a solution of the 
differential equations with the help of the method of 
difference . 

1 . Foppl 1 • Method 

The evaluation of the approximate solution for rec- 
tangular shell panels (fig. l) leads to the simple, jjrac- 
tical formulas below. Foppl 1 -s condition for the validity 
of the formulas, that the deflection of .the plate must be 
considerably greater than the sheet thickness is probably 
always complied with by the sheet thicknesses employed in 
airplane design. 

Deflection in panel center (point m) : 




stresses in panel center (point m)j 




stresses on the edge in the center of the short side b 
of the rectangle (point r-^) l, t 

"/I ? 
Cf z =n 4 / p E 

stresses on the long side of the rectangle, a (point r a ) 
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7 a a^ ... 
7r a W s 2 

v/here p is uniformly distributed hydrostatic pressure, 
kg/ cm 8 

£ modulus of elasticity, kg/cm 2 
a and b half rectangle side::, cm 

s sheet thickness, cm 
n x to n 7 coefficients 

The coefficients a, to n 7 were also determined 
for rectangular plates and can he read from figure 2 for 

any aspect ratio of the rectangle X = ~. 

a 

The maximum' stresses occur accordingly in the center 
of the edge of the long rectangle side in z direction, 
liven the square plate is stressed highest in edge center 
and not in plate center, as claimed by Foppl in his vol- 
ume "Drang und Zwang. 11 A new edition is to carry this 
correction. In order to convey a clear picture of the 
stress distribution over the whole plate, the stress coef- 
ficients n for the square plate are shown in figures 3 
and 4 plotted against the corresponding plate points. 

Example : 

Dimensions of plate, 60 X 60 cm 
Sheet thickness, s = 1.4 mm' 
E = 740,000 kg/cm 2 
Load, p = 20 t/m s 

Compute the deflection f in plate center and the 
maximum stress G on the edge. 



f = 0.8 a 3 /^ * 0.3X30X 7 §£22 = 1.99 cm 

v Bs v 740000X0.14 
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0.56 'JT(K-) 



' 2X30 



b 0.56 / 740000 ( ) = 3700 kg/cm" 3 

2. Hendry. 1 s Method 

This method was published in the Zeitschrift fur 
angewandte Mathemat ik und Meehanjk (1923.), pp . 81 and 
423 in the report entitled: "The Calculation of Thin 
Rectangular Plates with Vanishing Flexural Stiffness. 11 
For the square plate the difference method affords: 

/pa 

Deflection f = n i a / — 

v Es 

The ni factors can he taken from figure 5. 



tres ses : 3/ / \ * 3 



a = n o 
z 



The factors n 2 and n 3 are plotted in figures 6 
and 7, The determination of the factors for rectangular 
plates entails considerable raper work and has therefore 
not been made go far. 

For the same example as below 1 the deflection is 

0 7 P 

f = 1,99 — — = 1.8 cm 

0-80 . 



0 .436 

= **T* s 2700 ^7 = 2100 k s/cm 2 



C. COMPARISON OF EXPERIMENTAL RESULTS WITH 
THE CALCULATION METHOD 



For purposes of checking the cited calculation meth- 
ods as to practicability and reliability various 60x60 cm 
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plates of different thicknesses were stressed to failure 
under hydrostatic pressure. The employed test set-up is 
shown in figures 3 and 9. 

Tank and cover were of rugged construction in order 
to keep Interference of skin deflection in consequence of 
tank deformation to a minimum* 

The measurements disclosed a good agreement "between 
the computed stresses and deflections and the experimental 
results so long as pure elastic behavior of the material 
prevails. The principles of the theory are herewith 
realized. This stress range is decisive for our designs, 
"because according to the design specifications the strains 
under service loads (j = 1.0) must lie within the elastic 
limit and the permanent deformations therefore remain 
unimportant (5 percent). 

The experiments further proved pure elastic "behavior 
to "be tied to relatively narrow limits of the stress. The 
plastic form changes start very soon. In this range the 
calculation gives, of course, erroneous values, deflections 
too small, stresses too high. Sheets under tensile 
stresses, as we have here, arc capable of very pronounced 
plastic form changes and considerable permanent strains. 
Through this plast icizat ion of the sections the plate 
bulges out more, which tends to . effect a reduction in the 
stresses. The bearing strengtn can therefore be raised 
considerably above the mathematical values. The plastics 
theory lately has been concerned with the stresses in the 
plastic range, but the theory is still in its initial 
s tage . 



For confirmation of the foregoing arguments the test 
data of a 60x60 cm, 1.4 mm thick shell, plating are in- 
cluded, figure 10 shows the experimental deflections in 
the center of the shell panel plotted against the load in 
comparison with the mathematical values, along with the 
permanent deformations from the different loads after re- 
laxation to zero. It is seen that the plastic deforma- 
tions already start between 1 and 2 atmospheres and that 
in this area the theoretical and experimental deflection 
curves also disperse considerably. The plate does not fail 
until at 11.2 atmospheres; hence the elastic range up to 
2 atmospheres at most is quite small compared with the 
ultimate load. 
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According to calculation the following ultimate 
loads would have "been afforded for 0 2 = 4400 kg/cm 2 



Foppl : 




s 



2 / 4400 3 x0.14 i 



-j 



E a 2 v 0 . 56 3 x?40000x30 2 



= 3.8 atm 



Hencky 



4400 3 X0.14 2 ri . . 

p = / = = 7.14 atm 

" ,44° X74000L ^30 2 

The ultimate loads obtained are substantially higher 
for the stated reasons, 3P8ppl f s approximate solution com- 
pared to Hencky f s gives substantially lower values in the 
elastic range; hence it leaves one on the safe side. In 
figure 11 Hencky ! s computed elastic line in panel center 
is compared with the experimental for 2 atmospheres. 

The difference in this load stage is small because 
only minor permanent deformations are present, according 
to figure 10 (in the vicinity of the elastic range). 

The edge stresses recorded with tensiometer showed 
themselves too snail, since in consequence of the not al- 
together negligible flexural stiffness of the sheet the 
fixed end moment was still effective on the test station 
and to a lesser extent because the tensiometer recorded 
merely the length change of the chord but not that due to 
sheet curvature (fig. 12). The load recorded at 2 atmos- 
pheres in center of plate edge was 1700 kg/cm 2 , i.e., much 
lover than computed previously. 

When the stresses are computed from the measured de- 
formations, they are in better agreement with the calcu- 
lation. The elastic line (fig. 11) satisfies the parabola 
equation very satisfactorily 

x * if y(l - y) 
1 2 

so that the arc length can be ascertained according to the 
relation (fig. 13) 
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The total strain on length 1 is: 

= 4i = d - i a iff 

y l \ 3 M 

On the assumption that the stresses CT and cr r . 

over the axis of symmetry of the square are constant and 
that j z on the average is -jr 0" y (which is approximately 

justified on the 'basis of the n values of figures 6 
and 7) the elastic range follows Hooke's laws 

€ - i fa v - ^ - £5. 22. 

7 E V y m > 16 E 
V. r ith f bs 2.0 cm recorded in the present example, we have 



J y = if E e y = if 740000 I (f Y = 2700 k «/ 



cm 2 



The agreement with the computed stress according to 
FSppl is perhaps accidentally so good. The failure oc- 
curred in the center of the sheet edge, where the stresses 
in the elastic range also are maximum (fig. 14). 
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Figs. 1,11,12,13 




Fifmrs 1.- Notation of 
rectangular 
skin plating. 



Figure 11.- Computed and 
experimental 
elastic lines at 2 at. in 
skin panel center. 
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Figure 12.- Tensiometer 
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FigB. 2.3,4,5,6. 
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Figure 2.- n values for any aspect 
ratio - of rectangle. 
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Figure 5.- Hencky's values n^ 



Figure 6.- Stress factors n$ for ffg. 



J 



NACA Technical Memorandum No. 965 



Figs. 7,8,9,10,14. 
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Figure 7.- Stress factors n3 for Oy 
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Figure 8.- Diagrammatic sketch of 
test set-up. 
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Figure 14.- View of break of 

600x600x1.4 mm 
sheet at 11,2 at. 




Figure 9,- Testing device. 



Figure 10.- Deflection in panel 

center plotted against 
load, skin panel: 600x600x1.4 mm. 
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